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ABSTRACT 

This  work  adrssses  tha  application  of  digital  optimum 
control   theory    to   a    tank-to-turn    missiia. 

A  optimal  guidance  law  has  beei  developed  aad  tes-ed  in 
several  scenarios,  using  a  2-D  modsl.  Effects  cf  sample 
rata,  pitch  angle,  gravity  and  approximations  for  small  and 
larae    roll    excursions   are  discussed. 
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I-    INTRODUCTION 

Because  of  threats  from  highly  maneuvsrabls  high 
performance  aircrafts  and  the  need  for  increase  standoff 
ranges,  major  improvements  are  needed  in  guidance  and 
control   capabilities    of   missiles. 

The  high  maneuverability  of  targets,  has  Led  to  defense 
missiles  capable  of  develop  higher  Lift  accelerations  and  to 
more  complex  control  laws,  able  to  improve  performance  over 
well    know  laws    as   proportional   navegation. 

In  order  to  accomplish  these  nej  requirements  with  large 
standoff  ranges,  propulsion  systems  asing  airbreathing 
engines    has    been   studied   and   developed   in    recent    years. 

The  advent  of  airbreathing  engines  leads  natural  to  a 
consideration  of  bank-to-turn  missiles  in  order  to  minimize 
the    angle   of  attack    of  the    inlets. 

The  necessity  of  more  complex  control  laws,  leads  in  a 
general  way,  to  the  application  of  modern  control  and  -sti- 
mation  theory,  since  more  complete  informations  of  the 
states  of  missile  and  target  are  necessary  than  those  states 
informed  by  sensors  commonly  in  use  in  missiles  today.  This 
leads   to   the  use   of    a  airborne   computer. 

The  present  work  adr esses  the  design  and  evaluation  of  a 
optimal  digital  control  for  application  to  terminal  guidance 
in    a    fcank-tc-turn   missiles. 

One  continuous  two  dimensional  node!  was  adopted,  in  the 
following    form: 


X(t)=    A(t)    x(t)«-    B  (t)     u(t)+    E    g  (1.1) 
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where  the  effect  of  gravity  appears  explicitly  in  the  third 
term    on   the  righ   hand  side    of    expression    1.1. 

After  the  development  of  a  equivalent  discrete  model, 
the  optimal  control  problem  has  been  solved, using  a  modified 
Ricatti  equation  due  to  the  existence  of  the  third  term 
representing  to    the    gravity    effect. 

Next,  several  analysis  has  been  made  in  order  to  check 
the  effects  of  small  and  large  roll  excursions,  the  effect 
of  the  sample  rate  on  the  system,  and  the  effect  of  the 
initial  pitch  angle,  in  order  to  chsck  the  validity  of  such 
two  dimensional  model,  when  applied  in  some  scenarios  of 
inter  st. 
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II.    MODEL    OF    THE    SIS TEH 

A.  INTRODUCTION 

In  -he  present  work  the  problem  ot  terminal  guidance  for 
long  range,  bank-to-turn  missiles  »ith  ramjet  engines, using 
a  digitalized  system   has   been   investigated. 

The  model  developed  in  reference  1  is  used  as  the  base 
for  this  work.  After  the  digitaliza -ion  of  that  model,  an 
optimal   control    law    was   developed. 

B.  ASSUHPTIOMS 

Keeping   the    same    assumptions   as   in    ref.1,on=    has: 

The  missil  is  limited  :o  -2g's  and  +15  g*s  of  commanded 
acceleration  in  the  pitch  plane, with  zero  lag. Also  its  yaw 
auto  pilot  has  zero  lag,  yaw  regulator  maintains  zero 
sideslip. 

Missile  thrust    exactly    cancels  drag. 

The  angle  of  attack  is  assumed  to  be  v ery  small, which 
leads  to  the  commanded  acceleration  acting  normal  tc  the 
velocity  vector. 

The  missile  will  not  have  to  roll  through  a  large 
angle.  (Further  considerations  will  give  to  this  at  the  end 
of   the    derivation   of    the   control    law)  . 

C.  THE    CONTINUOUS    MODEL 

Using  the   same    reference   frames    as   in    ref.1,one    assumes: 
-Body   frame      with   xb    axis      parallel    tc      the   longitudinal 

missile    axis,         positive    y&    axis   out      of    the    left      wing    ,and 

positive    zb    axis    upward,  (see    fig. 2.1) 
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-Fligth  path  axis  with  xF  axis  parallel  to  the  velocity 
vector, positive  zp  axis  pointing  upwards  and  y..  axis 
pointing  to   the    left,  (see   fig.    2.2) 

In  fig  2.1  and  2.2r  the  angles  p  and  B  a:=  the  Eulerian 
roll    and   pitch    angles. 

The   state  vector    is   given   as 


x    =    [y...    yp    Aty     zp     zF  Atz     &ft  ] 


(2.1) 


where  y^  and  zF  are  the  :omponen-s  of  the  relative  target 
position,  yp  and  zp  are  the  relative  target  velocity, Aty  and 
Atz  are  the  components  of  target  acceleration , which  is  expo- 
nentially decaying    with   a   time  constant    <5  . 


A0  =  fl-k 


(2.2) 


where  ^   is   the    initial  roil    angle    (at    t=0) 
The    control    vector   is  given   as: 


u    =      Ac      Pc 


(2.3) 


where      Ac      is      the      commanded    acceleration      and      Pc      is      the 
commanded   roll    rate. 

The   nonlinear   plant   equation    is 


X     =f  (     X,     u    )      ♦ 


(2.4) 
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or 


X  = 


Aty  t  Ac  sin  d 
'Aty  /S 

Atz  -Ac  cos  6 

-au/g 

Pc 


0 
0 

o 
o 

-  J  cos  Q 

0 

0 


(2.5) 


whare   g   is      gravity's  acceleration   lad  9  is    -hs      pitch    a?.gl: 
as   seen    in    fig. 2.2 

Linearizing    and    set-ing 


G    =    E    q 


(2.6) 


on =    has 


X   - 


AtyhA!c(cosfio)Afi 
-Aty/<5 

z 

r 

Atz.  *  Ai(cosfo)&(p 

-Atz/ 6 

O 


0 

0 

smfc 

0 

0 

0 

0 

0 

-cosfc 

0 

0 

0 

0 

o 

(2.7) 


where 
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•-[ 


0      0      0       0      -ccs0        0      0 


] 


(2.3) 


in   aqn.    2. 7, we    hava    set 

cos^    =   cos(/4   +Aft)    =co3^    cos^-sin^    sir.  ^^ 

sin^zf    =   sin(/5    +  4^)    =sinj^    cos4/tc3  3^o    sin  4^ 

and    expanded  in   &fi    which   is   considered   small. 

Now      assuming      that      Ac, which    is      actually      the      desi 
control    Ac,    can    be    expressed    in   una    form    of: 


-«^ 


Ac    =    A  co 


Ti 


(2.9) 


with    Aco    =    Ac   at    t  =  0,      one    has 


x    =Ax   +    Bu    ♦   E    g 


(2.  10) 


wher  e 


A  = 


o      /      o      o 

OO/O 

o     o    -'/&    o 

OOOO 
O  0  0  o 
0  0  0  0 
OOOO 


O 

o 

0 

O 

0 

Aco 

'*■$ 

O 

0 

O 

/ 

0 

O 

0 

/ 

Aco 

>f 

0 

-*/3 

O 

o 


COS 


5in 


i 


(2.  11) 
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(2.  12) 


&    - 


0 

0 

j;o  f£0 

0 

0 

0 

0 

0 

-  CM  4 

0 

0 

0 

0 

1 

■  ■I 


0      0      3       0      -cos  a 


0      0 


(2.13) 


g    =    g 


(2.  14) 


D-       THE    DISCRETE    MODEL 

With  the  introduction  cf  i  ligital  ccmpu-er  to 
control  the  continuous-time  sys-em,:>ne  has  to  have  some  kind 
of  interface  in  order  to  take  cars  of  the  communication 
between  the  discrete  and  continuous-time  systems. In  this 
case  it  will  be  considered,  a--o-d  and  D--o-A  converters  as 
samplers   and  zero-order   holders    as   in    reference    2. 

In    such    case,   considering   -a=    system: 


x(t)     =    A(-)     x(z)     ♦    B(t)     u(t)     +    E(t)     g(t) 


(2.  15) 
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one    can    write  the  state    of    the    system   at   time   t(k  +  1)     as    : 

*<**>'    <\.,  »V  *^*J^\,  'TB(T  *|  U(V        <2-16> 

W/(tM     ,t|)     E("|)    d<|    g(tk) 


*« 


where  2$ (t, to)    is   the    transaction   matrix   of   the   system  repre- 
sented  by   eqn.    2.13. 

Furthsrmore,      we      will    considered   that      the   sampling 
instants   are  equally    spaced,    or: 


■       -  t     =  r  (2. 17) 


t         =    kT    ♦   I  (2.  18) 


so    one  can    replace 
t„     =    kT 


thus  , 


x(kT+T)=    ^(kT+D     x(kT)  (2.19) 


+     Azf(kT+T,/»|  )     3(i])     a<r|      jfkr) 
^  (kT*T,*j  )    E(<*|  )    d«j     g.(kT) 


KT 

or    in    a    simplified    notation: 
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x(k  +  1)  =  Ad  (k)     x(k)    ♦   Bd  (k)    u(k) 


(2.20) 


+  E4(M   g(M 


wh  3  r  9  , 


a<j  <k)  =  p  (kr*T#Jcr> 


(2.21) 


Ba< 


tfC  (kT+T 


7j  )    b  ('»i  )    a* 


k~ 


/k"l>T 
^(kT+r,«y|  )    B(*j  )    d^ 


kT 


(2.  22) 


(2.23) 


2.      Calculation   9.1  lh.±    Matrices    A  (k)  ,    3(k)    and   E  (k) 

It  is  straightforward  to  show,  using  the  sparseness 
of  the  matrix  A (t) ,  that  ths  transicticn  matrix  of  equation 
2-19    is    : 


$(*t+t)  = 


/ 

T 

Ad 

A  5 

0 

0 

0 

M,, 

0 

f 

H.S 

0 

0 

0 

^z.r 

0 

0 

eT/* 

0 

0 

0 

0 

0 

0 

0 

f 

r 

H, 

Hr 

0 

0 

0 

0 

/ 

H< 

K 

0 

0 

0 

0 

0 

e* 

0 

0 

0 

0 

0 

0 

0 

/ 

(2.  24) 
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where: 


Ad       =Ad      =  5    T    '  5  2 

l.i  4.4 


.TV 

(1-e  '« 


) 


(1-e"%  ) 


For   the    calculation    of    the   others   terms    one    may    mak< 
use    of   the    property    of  the    transiction   matrix   that: 


dt, 


=    A(t    )    Ad(t     ,t(  ) 


so, 


dAd(kT+T(  |cT) 

=    A(kT+T)     Ad  (kT+T,kn 

dCkTi-T) 


For   A    (2,7)  : 


d  Ad„(kT+T,kT) 


=   AJ7  (kl+T)     =    Aco 


,_  A. 


~os£ 


d(kTtT) 

fL&Lw.CG:L.!£tL«  aco 

d(kT*.T) 

/kT»T 
["  1-  -^-r--^-  J     a(kT+T) 

H.,=  Aco  «34   [T-[-yrlJ-]    I2  j 

For    A    (1 ,7)      : 
dAd,.r(kT+Ty  «t) 


"•ri 


kT+T 


d(kT*T) 


ld2j{kT+T!,  kT) 


ieT*T 


Ad     (kT+T,kT)     =         Ad       d(kT+T) 

KT 
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&d-,r  [T2-[-W~J  T3]  A"  cos^ 


Doing   the   same  process    for   Ad  (5,7)       and   Ad  (4,7)       one 


has: 


=    Aco   s: 


For   th  e      derivation    of   the      matrix    Bd  (k)         one    needs 
according  to  egn.    2-22 

ft  (kT+T)     B(1|)=    Ad     (kT+r,»|)     B  (  1|  ) 

where, 

/    A*     A*     0     O     O    A* 
0    /      An      0     O      0    A« 

'2,  j  '2.r 

(9      0      At,  0       0       0        0 

0  0  0  /  Am  Am  A* 

0  0  0  0  t  Am  Am 

(9  0  0  0  0  4~  0 

(9  0  0  0      0  0  i 


A<j(kT+T,*]  ) 


where    A*  represents    Aj(k?+r,-ri  ),    and 


fi^J   = 


0 

0 

sin  d 

0 

0 

0 

0 

0 

_  cos  £> 

0 

0 

0 

(9 

0 
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thus 


AdOcT*T,<rj)ftfr|)  = 


A\2s,ni 

A\r 

5.  a  $0 

A«U, 

0 

0 

-A\r*k 

AV.i 

-cos    fc 

0 

Am 

0 

0 

0 

whsre 


l1u"    kT  +  T"  1    =  *"U, 


A  ir\ 

Aro 
'»,7 


[* 


A CO     CCS 


(kr  +  T)    - 


/Ji.tJ  — 

~Yff" 


(kr+r)    r2  - 


-    T  T)     * 


T2 


1 J 


Ace    cos 


'\r 


T(kH-T)     - 


JJl+J. 
"TtT  " 


(kH-T)     T2    - 


-   T  Of)    ♦ 


kT+r- 


JJltl. 

(kTfT) 


T2 


AZD 


ln/l 


'5.7 


Ti  '  ^T.-       J  r° 


and    for    B ,  (k) : 


22 


b. 


B 


C»M 


&d 


im  po 

5d 

0 

0 

H, 

Ba 

-  Co5    tf£ 

5d 

0 

0 

£ 

T 

1.2 


hi 


4.2 


J.* 


(2.25) 


where 


Bd 


,kTtr 

,,=  /    (**♦*"  «|)     d1     sin£    =   -I-   3in/o 

KT 


T' 


/  kT  +  T 


which    after    some    algebric  work    has    been    found 


3        (k)  = 


Z 


(kT+T)  (kT+t)     kT  (kt) 


l    1 


3T\ 


Z  Ti 


6Ti       J 


ace   cos  (6 


and 


kT+T 


BVk)  =  /    %a1 


kT 


which    can   be  found    to  be 
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Bd1<t(M  = 


2T\ 


"J4 


Aco    COS 


thus 


Bd-,0O  = 


*dMW- 


1 

7 


T* 


A  CO     S 


i»/. 


3  Ti  2  T;  6T) 


J  -, 


A  co   siniD 


In    the   same    way: 


0(kT+T,y)    E(<r\)  - 


0 

0 
0 

-t co6  e 

-  cos    0 
0 

o 


thus    E    (k)     is   squal    to 


£       SB 


0      0      0 


T4 


cos  6      -I    ros  6) 

2 


0 


where      we      have    considered   B   as      a      constant      angle.  (Further 
comments   en   this    after   development    of    the   control    law). 

Notice  that  througout  this  work  rthe  commanded 
acceleration  has  been  considered  an  unknown  and  the  assump- 
tion has  been  made  that  it  will  be  in  the  fora  of  eqn.2.9. 
One  will  need  in  further  developments  to  consider  the 
control  Ac  as  a  known  Ac  (k)  ,  which  will  be  a  constant  between 
kT  and  kT+T.  With  such  assumptions,  the  discrete  represen- 
tation  of  the   system    is    easily    found    to    be    : 
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x(k  +  1)  =   A  (k)    x  (k)  ♦    E(k)    a  (k)  ♦   E   g 


(2.26) 


whsre 


*d    = 


/ 

0 


0 

0 


T      /lei      0      0       0      A<* 


/.j 


',7 


/  Ajiit    0      0      0     Adzjr 

0  <fT/s    0^00 

0      0  0/r     ^   H,, 

0      0  0       0/     Ad^    A^r 


0      0       0       0     c 


r/s 


0      0       0       0     0 


0 
/ 


(2.27) 


with 


-T 


Ad    =    3T-32     (1-    e  '6  )    =    Ad 


r/G 


•»ll 


4,6 


■T/s 


Ad     =      3(1"    9    /6  )     =    Ads. 

AdliT=    "2"    AC    C0SPo 
hd     -   T    Ac    cos/o 
Ad4>7=   .11    Ac   sin/0 


Ad57=    T    AC    sin/o 
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and 


Bd(K) 


7-    S,n/o 

-1-   Ac   cos  p 

T    si n  ^, 

y     Ac  cos^ 

0 

0 

-     X.     COS  C6 

T3                        / 
—    Ac   5inf6o 

-  T  cos  0 

'  0 

—      4c  sin  d 

0 

0 

0 

r 

and 


Ed(k) 


0 

0 
0 

_  Ji  cos  e 

z 

-   T    cos    0 
0 
0 
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Zf    1 

i 

Y 

b 

Ac  C0S(?3 

r         "71 
0/ 

/  Zb 

'Ac 

1 

f 

Ac  stn  <fi 

Figure  2.1    Reference  Frames. 
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4  Z 


V=Xb 


Figure  2.2    Referents  Frames. 
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III.  THE  OPTIMAL  COHT ROLLER 

A.   DERIVATIOH  OF  THE  OPTIMAL  CONTROLLER 

In  order  to  have  a  suitable  guidance  law  ;o  implement 
the  control  commands, we  will  liaimizing  the  following 
performance  index: 


H.I 


J=  -j-  xr(n)H  (n)x(n)  +  )  ---  ur(k)  2  <M  u  (k) 


(3.1) 


*»o 


where   x  (N)  is  the   final  state   at   t  =  ri 

As      we    want      to      minimize      the   final      miss      distance, the 
weighting  matrix    W(N)    is   taken   as 


/        O       0       0       0       0       o 


0       0       0       0       0       0 


ooooooo 


0  0  0  f  o  o  0 
ooooooo 
ooooooo 
o       0      0      O       0      o      o 


(3.2) 


and      <2(k)         is    a      two     by      two   positive      definite      symmetric; 
weighting   matrix   to    be  chosen. 

In   the    derivation  of   the   solution.,      reference    3    has   been 
followed   keeping      in    mind      that   the      state    sguation      has   the 
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x(k+1)=  A  (k|  x  (k)  +  B(k)  u(k)  +  E  j 


(3.3) 


or 


x(k  +  1)  =  f  (x(k),u(k)  rg) 


(3.U) 


where   "he  -third   term, which  represents   -he   effect  of   the 
gravity  has  beer,  considered  constant. 

Considering  that  the  performance  index  is  in  the  form: 


M-l 

J  =  ^[x(N)  J  ♦  \   L(k)  j"  x(k)  ,u(k)  ,g(k)  J 


(3.5) 


k  =  0 


we   need    to    find    a   sequence    of    u(k)     that    minimizes    J. 

Adjoin    the    system  equation   to   J    with   a    multiplier    A  (n) 


At-' 


=  ^[x(N)j 

*    A  (k+1)     J  f, 


x  (k)  ru  (k)  ,g 


-    x(k  +  1) 


) 


and    defining  a    scalar   sequence    H  (k) 


H(k)=    L(k) 


x(k)  ,u  (k)  ,g 


+  A (k+1)     fK  j"   x(k)  ru(k)  ,g] 


k=0,  1 ,2, n-1 


(3.6) 


(3.7) 
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one    has: 


J=  f6\  x(N)l    -  A(N)     x(N) 


N-l 


*      \    [h  (k)    -  aV)     x(k)    J    +    H(0) 


Ks  I 


Considering    differencial   changes    in    J 


(3.8) 


dJ  = 


-£-   -     A   (N)         dx(N)   + 


L    £*(») 


(3.9) 


N-l 


«^u(k)  ^x(o) 


♦  ^&    du(0) 

^U(O) 

choosing  the  nultiplier  A(k)    so   that 


c^x(lO 


(3.  10) 


thus 


^x(k)  a>X0e) 
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and 


o>H_0O 

Jo(k) 


=    0 


(3.  11) 


with    boundary   condition 


X(N) 


(3.12) 


we   obtain  the   minimization    of   the   per formance   index. 
In   the   present   case    we    have: 


«-i 


J=    y   xT(N)    »(N)    x  (N)    ♦  \     uT(k)    2(k)     u(k) 
+     A  (k+1)   I  A  (k)  x  (k)+3(k)     u(k)  +E  3-x(k+1)l 


(3.13) 


and    H  (k)  : 


H(k)    =   y    uT(k)    a  (k)    u(k)  ♦ 


(3.  14) 


♦    X   (k+1)       A(k)    x(k>     +    3(k)    d(k>     +    S   g 


then    in    order   to    minimize   J: 


----   =   u    (k)     Q  (k)     ♦    X  (itH)     3(k)    =0 

o)o(tf) 


(3.  15) 
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or   considering    that    Q  (k)  =3  (  k) 


Q  (k)    u(k)    =  -    B    (k)    >    (k+  1) 


and 


o)x(k) 


(3.16) 


A(k)    =        )J"(k+1)     A(k) 


and 


X(N)     =    x'(N)     H  (N) 


(3.  17) 


Notice      that    we      are  not    weighting      the      states    in      the 

psrfcrmance   index,    except  the    last   state. A    mors    general  form 

could   be   obtained,         with  ill    s-atas    baing      weigthing    ,  if   wr 
change   -he    sqn    3.16    to: 


\(k)    =        >(k+1)     A(k)    ♦    x    (k)     W1(k) 


(3.13) 


whars    W1  (k)    is    the    weighting   matrix    of   the   s-r.as. 

With  equations  3.15,3.16  and  3.17  ona  is  abla  to  find 
tha    sequence   t>f    u(k)     that    will    giva    th3    minima   controls. 

Such  set  of  aquatics  can  be  solved  by  the  sweep  method 
as    in   ref.2    . 

We   will   look   for    a  solution   of   tha   form: 


u(k)=-F(k)x(k)  -FG  (k)g(k) 


(3.19) 
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what    msans     that    the    commanded      acceleration   and      roll    rate, 
will    have  a   correcti.cn  due    to   the    effect   of   gravity. 
Placing: 


A(k)    =    S(k)     x(k)     +  SG(k)    g  (k) 


(3.20) 


from    egn.    3.15 

Q  (k)u(k)  =-B    (k)  |S(k+1)  x(k-M)  +53  (k+1)  g(k+1)J 

from    egn.    3.3 


3(k)    u(k)    =  -B    (k)     S(k+1)       A  (k)     x(!c)  t 


[ 


(3.21) 


♦B(k)    a(k)    ♦   S    gj   -   BT  (<)    SG(k+1)    g(k  +  1) 


so 


]■ 


Q  (k)    +   B    (k)     S(k+1)    3  (k)       u(k)  = 


-BT  (k)     S(k  +  1)     A(k)     x(k)    -    B     (k)     S(k  +  1)     E    g  (k) 


-BT  (k)     SG  (k+1)     g  (k  +  1) 


considering  that    g    is  a   constant 


u(k)    =-[Q(k)     +    BT(k)    S(k+1)    B    (k)J 


(3.22) 
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B     (k)     S  (k+1)    A  (k)     x(K)  + 


[BT(fc 


)     S(k  + 1)      E  +    B    (k)    Sg  (k+1) 


] .] 


so: 


whsre 


u  (k)  =  -F  (k)  x(k)  -FG(k)  g 


■[ 


.1 


F    =      Q(k)     ♦   B     (k)     S(k+1)     B(k) 


(3.23) 


.    B    (k)    S  (k+1)     A(k) 


FG    =   |q(M    +    B    (k)     S(k+1)     B(k) 


'] 


(3.24) 


[3T(k)    S(k  +  1)     E  +    BTSg(k  +  1)l 
from    eqn.3.16    and  3.19 

X(k)=AT(k)     A(k+1)=AT(k)    |  S(k+1)  x(k+1)  +  SG  (k  +  1)  g  (k+1)j  z 
=AT  (k)     S(k+1)    A(k)     x(k)     +    AT  (k)     S(k+1)    B(k)     u(k) 
+  AT  (k)     S(k+1)    2g    +    AT(k)     SG(k+1)     g(k+1) 
from    eqn.    3.19 

X(k)    =    AT(k)    S(k  +  1)     A(k)    x  (k)     +    A*  (k)    S  (k  +  1)     B  (k)  . 
•[-F(k)     x(k)     -   FG(k)     gj    +    AT  ( k)     S(k+1)     Eg    + 
+  A    S3 (k)    g  (k+1) 
so,    as   g   is    a   constant: 

S  (k)x  (k)  +SG(k)  g=      AT  (k)  S  (k+1)  A  (k)  -A  (k)  S(k+1)B  (k)  F(k)~|  + 
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+\aT  (k)     S(k*1)    E   -    AT(k+1)    3<k)     FG(k) 


♦    A    (k)     SG  (k+1)  J  g 

:hus 


S(k)=A     (k)S  (k*1)  A  (k)-A    (k)  S(k+1)  B(k)  F  (k)  (3.25) 


and 


SG  (k)=AT  (k)S  (k  +  1)  E-AT  (k)  S  (k+1)  B  (k|  FG(k)  +  (3.  26) 

♦    AT  (k)     SG(k+1) 

These   equations,    3.25f     3.26,    3.23    and   3.24   can    be    solved 
backwards    with    the    final   cor.di-ion: 
S  (N)  =W(N) 

SG  (N)  =0 

Notice  that  this  satisfies  our  previous  boundary  condi- 
tion   in   eqn.    3.17   where: 

(N)  =  W(N)  x  (N) 

S  (N)  x  (N)  =  W  (N)  x  (N) 

so 

S  (N)=W(N) 

B.        EXTENTIOM    OF    THE    MODEL    FOR    LARGE    ROLL    ANGLES 

Up  -o  this  point,  one  has  to  take  into  account  tha- 
throughout  the  development  of  this  work, the  angle  has  been 
considered   small,    it    is   necessary   ts    ralax   this   rss-riction. 
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In  order  to  do  that,  the  system  has  been  broken  in  two 
blocks   as   in   figure    3.1. 

The  first  block  is  a  representation  of  the  algorithm 
which  will  calculate  the  optimal  commands.  The  algorithm  has 
contain  with  itself  an  exact  model  of  the  system  or  missile. 
The  model  of  the  missile  is  initialized  from  the  information 
on  the  initial  states,  the  initial  input  command  Aco  and 
initial  roll  angle  (at  t=0  ) ;  and  computes  the  optimal  gains 
and  further  the  optimal  commands  which  will  be  feed  to  the 
missile. 

The  method  adopted  in  computing  the  optimal  commands  is 
more    easily   understood  if   one   considers    figure   3.2. 

In  figure  3. 2,  the  lines  numbered  as  0  in  the  graph  for 
Ac  and  for  Pc  are  the  optimal  commands  for  a  given  initial 
roll  angle (0Q)  .  Lines  number  1  are  the  commands  for  a  second 
initial  roll  angle(j^)  larger  than^,  and  so  on.  Thus  in 
figure  3.2  one  has  a  family  of  optimal  commands  for  any 
initial   roll  angle. 

Notice  that  the  upper  line  of  the  graph  of  Ac  represents 
the  accelerations  of  a  missile  whioh  had  at  t=0  a  correct 
initial  roll  angle  in  order  to  hit  the  target  with  no 
comands   in    roll    rate. 

In  the  present  method  the  computer  performs  the  calcula- 
tion of  the  commands  only  for  the  first  step  of  time  and 
then  feeds  these  commands  to  the  missile.  The  missile  is 
then  driven  to  the  next  state  (x(k+1|)  and  feed-backs  to  the 
computer  the  information  on  the  roil  angle  at  that  step. The 
roll  angle  feed-backs  from  the  missile  is  considered  by  the 
algoritm  as  the  initial  roil  angle  at  t  =  0  and  the  next 
commands  are  calculated. Not  ice  that  at  this  second  step  the 
algorithm  will  feed  to  the  system  the  second  command  (at 
t=t,  )  .  This  process  is  them  repeated  until  t  is  egual  to 
the    intercepted    time. 
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It  is  important  to  realize  that  with  -his  method  of 
calculation,  since  the  algorithm  was  developed  with  the 
assumption  of  small  roll  excursions  some  error  is  expected 
due  to  the  fact  that  in  oomputing  the  gains  by  solving  a 
Ricatti  equation  backwards,  as  has  been  done, it  is  necessary 
to  update  the  system  from  t=0  to  t=ri  at  each  step,  and  in 
this  process  the  roil  angle  is  not  small. Notice  however  that 
we  are  applying  the  commands  only  in  one  step, and  if  one 
expects  that  the  roll  rate  will  decrease  to  zero, as  we  are 
increasing  in  time,  the  variation  of  the  roll  ar.gls  will 
tsnd  to  decrease,  so , we  can  expected  that  the  error  will 
decrease   as   the    time    increases. 

Another  important  point  to  be  studied  is  how  the  missile 
itself  (second  block  in  fig. 3. 2  )  has  to  be  implemented  in 
order   to   be   valid  for  small    and   large   roll    excursions. 

It  is  considered  that  one  has  the  perfect  knowledge  of 
the  commands,  thus  the  missile  is  modeled  as  a  state  vari- 
able system  as  in  eqn  2.26  with  the  initial  roil  angle  being 
update  at  each  step.  In  this  way  the  system  will  take  the 
initial  roll  angle  as  the  summation  of  all  previous  initial 
roll    angles,    (see  fig  3.4) 

From  the  original  variables  one  has  for  the  state  &6  , 
the    following: 


A^k*1)=A^(k)«-    ?c(k)     T  (3.27) 


which    is   shew    in   figg.      3.3,    where    the    initial   roil    angle   is 
kept    constant   and    Ap  is   update  each    step. 

Notice   however      that   considering      large   roil      excursions 
(see    fig.       3.4)       and    keeping   in      mind    that    the   angle    p      has 


been    defined  as 


38 


the    expression      3.27    is      not   valid,         since    in      modeling  the 
system   it   was   assumed  that     &p  would    be    samll. 

This  leads  to  a  change  in  the  expression  for  the  varia- 
tion of  A$  as  in  fig.  3.4.  In  fig.  3.4,  the  angle  &</>  is 
update    at   each   step,    so   one    has: 


^(k+1)  =   Pc  (k)    T 


(3.28) 


By  consideration  of  the  equation  2.27  it  can  be  seen 
that  by  setting  the  element  A(7r7)  to  zero,  one  can  obtaim 
equation    3.28. 

Thus  the  missile  model  for  large  roll  excursions  will  be 
represented   by    tha    following   equation. 

x(k  +  1)  =   A  (k)     x  (k)  ♦    B(k)    u  (k)  ♦    2    g 

whsre 


Alk)  = 


/ 

T 

\> 

0 

0 

0 

\, 

0 

1 

-4*1 

0 

o 

o 

*., 

0 

0 

-r/c 
e 

0 

0 

o 

0 

0 

0 

0 

J 

r 

A«.i. 

A>., 

0 

0 

0 

0 

/ 

A,.. 

\, 

0 

0 

0 

0 

0 

f* 

0 

0 

0 

0 

c 

0 

0 

Ar< 

(3.29) 
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where    A(7,7)=0    ,for    large  roll   excursions 
A(7,7)  =  1    rfor    small   roll   excursions 


B 


"J    S>«tfo 

—    Ac    cos  d> 

T     Sif\   & 

—     Ac    cospa 

0 

0 

2               ^° 

—     Ac      Sin^o 

•  T   cos    & 

-I-      Ac    5">^o 

0 

0 

0 

T 

(3. 30) 


0 

0 

0 
II  cos  e 

z 

r   cos  0 

0 
0 


(3.31) 


We  will  redefine  the  states  x,  ,  as  zha  relative  posi- 
tion, x2  as  -he  relative  velocity  in  the  Y  direction  and  x3 
as  the  -arget  acceleration  in  Y  direc-ion.  The  stares  x^  , 
xs   and   xfc    has   the  same   meaning,        bat   in   the    Z   direction    and 
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xT   is  &p  .        In    this   representation   of    the   system,    x,  (k)      i« 
equal    to      the   component      of    the   miss      distance  along      the   Y 
direction   and: 


x,    (k+1)=    yp(k*1)  =    y„(k)  +    yp  (k)     T+  (3.32) 


♦  fc  I-  52  (1-  e  ^  M  Aty(k)  ♦-"£-    Ac(k)  cos^    4)Z*(k) 

♦  -j-    siniO   Ac(k)  +    ---    Ac(ki     cosp       Pc(k) 

The  first  three  terms  in  the  RH3  of  ear..  3.32,  are  easily 
seen  as  the  contribution  to  the  miss  distance  of  respec- 
tively the  previous  miss  distance,  relative  velocity  and 
target  acceleration.  The  following  two  terms  represents  the 
contribuiticn  cf  the  commanded  acceleration  and  the  last 
term  represents  the  effect  of  coupled  Ac  and  Pc,  and  tends 
to    be   small    due    to    the  cube    cf    the    sample    period. 

For  the  component  of  miss  distance  in  the  Z  direction, 
one    has: 

i    (k+1)«    zp(k4-1)=    z_(k)  +    zF  (k)     T  + 

♦  [5  T-  5  2  (1-  e"%  )       Atz  (k)  +    -j-    Ac(k)    sin^^^(k)    - 
-    ---    cos  !*£  Ac(k) +•    ---    Ac(k|     sin^     Pc  (k)  -    -j-    ccs  9  g 


where  its  ■'■erms  have  the  same  physical  meaning  as  in  the 
expression  for  x  (k+1),  with  the  effect  of  the  gravity  added 
to   the    expression. 

Since  one  can  notice  that  in  tae  representation  of  the 
miss  distance  in  T  direction  appear  two  terms  as  a  function 
cf  cos  ^  ,  and  in  the  representaction  of  the  miss  distance 
in    Z    direction    appear  two  terms   as   function    of  sin^    ,    it   is 
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interesting  to  verify  that  the  fourta  term  in  the  RHS  of 
both  expressions  acts  like  a  correotion  for  the  fifty  term. 
Befering  to  fig.  3.4b,  one  can  see  that  at  any  step  of  rime, 
the   commanded   acceleration    is   actually, 


Ac   ccs po  -  A  Ac 
and    considering    small  angles: 
A  Ac=   Ac   cos ]60  -   Ac   cos  ( j60  +  Afi)    - 

-   Ac  cosj60    -   Ac     <zos p  Ap-   sin  fa  sin  4p/  = 
=   Ac   sin^zf    A$ 

The  same  idea  can  be  applied  to  the  expression  for 
x    (k  +  1)  . 

The  terms  x^  and  xs  ,  represent  the  relative  velocity, 
and    are: 

-V 

x2  (k+1)=    y    (k)    +    5    (1-    e  /c  )     Aty(fc)     *    Ac(k)     cos^     T 

♦    Ac(k)    sin^     T   +    ---    Ac(k)     cos^    ?c(k) 

x5(k*1)  =    zp(k)    +    5    (1-   e"/c  )    Atz  (k)     ♦    Ac(k)    sin^    T    - 

-    Ac(k)    ccs^o   T   ♦    -y~    Ac  ( k)     sin^o    ?c(k)     -    T   cos  B     g 

Where  the  two  first  terms  in  the  RHS  represents  the  effect 
of  the  velocity  and  acceleration  at  a  previous  step,  and  the 
other  three  terms  has  the  same  meaning  as  previously  stated. 
The  terms  xa  and  x6  are  the  target  accelerations,  in 
this    model    being  exponentially   decaying. 


1  •      Effects    on    the   Mis s    Disiaaca    of    the    Extent ion    of   the 
Model 

In  previous  subsection,  a  extantion  of  the  model  for 
large  roll  excurssions  has  been  performed.  Notice  that  there 
are    two   models    of  the  system   being   used.    The    first    one,    used 
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in  the  algorithm  is  valid  only  for  samll  roll  excursions, 
and  a  second  model,  valid  for  small  and  large  roll  excur- 
sions  used   as  a    representation   of   the    missile. 

The  algorithm  with  the  first  model,  as  explained 
before,  is  initialized  at  each  step  with  the  actual  roil 
angle  of  the  missil,  and  performs  the  calculation  of  the 
comands. 

In  order  to  check  the  effect  of  the  extention  of  the 
model  on  the  miss  distance,  one  cm  define  a  ideal  initial 
roll  angle  (  ^0(oai  )  ,  as  the  roll  angle  at  t=0  in  order  to 
have  the  comanded  acceleration  vector  pointing  to  the 
projected  final  target's  position.  This  means  -hat  the 
missile  would  not  have  to  roll  tc  hit  the  target  (see 
fig3.5  ),  thus  the  commanded  roll  rate  calculated  by  the 
algorithm  will  be  equal  to  zero.  This  implies  that  from  that 
point  ahead,  the  roll  angle  is  constant  and  equal 
to  0O  ltita(  ,  and  that  the  time  history  of  the  control  Ac  will 
be   a    straight    line. 

The  fact  that  the  roll  angls  will  tend  to  this  limit 
deserves  an  investigation.  Notice  that,  as  show  in  fig.  3.2, 
at  the  moment  the  missile  reachs  its  maximum  roll  angle,  the 
control  Ac  will  be  the  required  acceleration  to  hit  the 
target   if   the   missile  initial   roll   angle   was    d0  .        This 

means  that  the  control  Ac  computed  would  be  correct  only  if 
the    missile   would   have  turned    immediatly    tc    this    angle. 

For   very    small  angles,    the   previous    comment    would   be 


acceptable,        but      in     a   normal      situation, 


showed 


figure3.5,  the  missile  will  only  reach  the  ideal  initial 
roll  angle  after  some  time,  and  due  to  the  vertical  compo- 
nent of  the  acceleration,  when  this  occurs  the  missile  would 
be  in  a  position  abcve  the  ideal  trajectory,  which  ireans 
that  it  would  follow  a  course  parallel  to  the  ideal  trajec- 
tory   to    intercept. 
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One  can  see  that  such  problem  will  lead  to  a  large 
miss  distance,  in  the  case  that  the  target  acceleration  it 
is  not  small.  Thus,  some  correction  is  necessary  in  order  to 
improve   the   missile    performance. 

Figure  3.5b  shows  the  missil  at  some  point  of  its 
trajectory  where  it  has  reach  its  maximum  roll  angle,  thus 
it  is  at  a  parallel  course  with  its  ideal  trajectory  to 
intercept.  it  this  point,  since  ill  the  states  are  known, 
it  is  possible  to  recompute  a  new  <rf idea\  •  So,  if  at  this 
point  the  computer  is  feed  with  with  the  states  at  this 
point  it  will  compute  the  commands  in  order  to  drive  the 
missil  to  the  new  $taeo\  »  which  will  introduce  the  desired 
correction.  Notice  that  such  correction  can  be  made  during 
all    the    flight,    from    t=0 ,    to  t=Ti. 

In  the  present  work  this  method  has  been  accom- 
plished by  feeding-back  to  the  computer  the  roil  angle  and 
all  the  states  of  the  missile.  With  this  information  the 
computer  is  able  to  perform  the  calculation  of  the  corrected 
commands  at  sach  step  of  time.  howeveras  the  states  are 
being  updated,  it  is  also  necessary  to  update  the  time  to 
intercept,  which  has  been  done  using  the  time  to  go  to 
intercept,    or: 

Ti  (k)  =   Ti-  k    T 

where  Ti   is   the    nominal   time   to   intercept. 

C-        SIHULATION 

In  order  to  keep  the  same  assumptions  as  reference  1, 
the  matrix  Q  (weighting  the  control)  has  been  put  as 
suggested   in  such   reference    or: 
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b, 


(3.33) 


with 


b.  =5.78      103 


b2=5.0   meters 


Five   diffarent    cases   ware    run: 

Case  0,  tested  with  one  simple  medal  valid  only  for 
small  roll  angles  begins  with  missil  and  target  on  parallel 
courses  to  tha  inert ial  z  axis,  tha  target  100  meters  abova 
the  missile  and  with  an  evasiva  manoavar  exponentially 
decaying  with  time  constant  of  20  seconds.  The  initial 
acceleration  of  target  was  -.5  g*s  in  y  direction,  (sea  fig. 
3.6) 

Case  1 ,  with  the  same  scenarios  as  case  0,  but  was  run 
using   the  algorithm    for   large   roll   excursions. 

Case  2 , is  the  same  scenario  as  in  case  1,  except  that 
tha  initial  acceleration  of  target  was  -1.0  g's  in  •'■-he  y 
direction.  (Same    as   case    1    in   raf .    1|  . 

Case  3,  the  same  as  case2,with  target  acceleration  of  -4 
g's.  (   Same    as   case    2    in   tha    reference    1) 

Case  4,same  as  case  3  but  with  target  at  initial  posi- 
tion 600  meters  bellow  the  missile. (S ame  as  case  3  in  refer- 
ence   1.) 

D.       COMMENTS    AND    CONCLUSIONS 

1  .       Results 

In  case  0,  the  missile  begins  its  trajectory 
commanding  26.5  m/sec2  and  the  time  history  of  the  control 
Ac    follows    exactly    a    straight    line      in    the    form    suggested    in 


45 


ref.  1 ,  as  seen  in  fig. 3.  7The  control  Pc,  baginir.g  at  .35 
rad  per  second,  is  decaying  and  raachs  zero  at  k=80  (see 
fig. 3.8).  Figures  3.9  and  3.10  shew  the  miss  distance,  where 
one  can  see  that  the  missile  is  crossing  the  target  with  a 
CG-to-CG   distance   of    1.5   meters. 

Figurs3.11,  shows  the  tiae  history  for  the  roll 
angle,  which  as  expect,  reachs  a  :Dnstant  value,  with  the 
missile   crossing   target    at    t=Ti,    with   a    bank    of    .44    radians. 

In  case  1, which  was  run  with  the  model  for  large 
roll  excursions,  the  missile  kept  the  same  Aco,  but  there  is 
a  very  small  increase  in  further  commanded  acceleration  in 
order  to  correct  the  effect  of  the  roll  angle  on  the 
vertical  component  of  the  control  Ac,  as  seen  in  fig. 3.  12, 
and    table   I. 

Figure  3.13  shows  that  tie  roil  rate  decreases 
almost  as  before,  and  the  final  roll  angle  is  .42  rds,  as 
shown  in  fig. 3. 15  and  table  I.  Refecing  to  figure  3.14  there 
is  a  change  in  the  final  miss  distance,  which  is  better  than 
case  0,  due  to  a  improvement  in  its  3  component  (see  -able 
I).  This  results  in  a  final  CG-to-C3  miss  distance  of  .65 
meters. 

In  case  2,  the  missila  has  tha  same  Ac  at  -=0,  but 
with  the  corr action  for  roll  angls  being  increased  due  to 
the  increase  of  target  acceleration  (see  fig. 3.  16),  the 
commanded    Ac  raachs    a   larger   peak   value. 

The  initial  value  of  the  commanded  roll  rate  is  .69 
radians,  which  is  larger  than  case  1,  due  to  the  increase  in 
the      target   acceleration.  In    figure      3.13,         -here    is     nc 

noticeable  change  in  the  shape  of  tha  ourve  for  Z  direction, 
and  in  the  Y  iir action  the  final  distance  is  about  the  same 
as  in  case  1.  The  CG-to-C3  distanoa  at  t=Ti  is  .73  meters. 
The  larger  roll  rate  leads  to  largar  roll  angles  as  seen  in 
fig.     3.19,    where   the    final    roll   angLs    is    .72    radians. 
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The  effect,  of  target  acceleration  can  be  easily  seen 
in  case  3,  where  one  can  see  that  with  the  same  Aco,  the 
accelerations  are  largelly  increased  from  this  point,  and 
the  missil  begins  its  trajectory  with  very  high  roll 
rate (see  figures  3-20  and3. 21) .  There  is  a  change  in  the  Z 
component  of  the  miss  distance,  that  decreases  its  final 
value  to  .05  meters,  but  now  the  aiss  distance  in  Y  direc- 
tion is  made  worse  as  shown  in  figure  3.22,  which  leads  to  a 
final  C3-to-C3  distance  of  1.5  metsrs.  In  fig. 3. 23  one  can 
see  that  the  missile  crosses  the  target  with  a  bank  angle  of 
1.26    radians. 

In  case  4,  due  to  the  position  of  target  600  meters 
under  the  missile,  the  initial  commanded  acceleration  is 
negative  and  reachs  the  limit  of  -2  g's.  The  initial  roll 
rate  begins  at  a  smaller  value  than  in  case  3  but  increases 
during  the  initial  part  of  the  flight  reaching  its  peak 
value  at  1.75  seconds  when  again  as  in  the  previous  cases 
begins  to  decay.  As  the  missile  banks  tc  roil  angles  larger 
than  90  degrees,  the  acceleration  goes  to  positive  values, 
as    seen      in    figures      3.24   and    3.25      Figure    3.26,  shows    the 

worse  case  among  these  in  respect  to  the  miss  distance, 
mainly  in  the  Z  direction,  and  in  the  final  cg-to-ca 
distance,      which   is      equal    to    4.43   neters.  Also   the    final 

roll  angle  of  3.0  radians  is  the  Largest  among  all  these 
cases,    as   seen    in   fig.    3.27. 

2 .      Comment  s 

Defining      the   projected      zero      effort   miss      distance 

(ZSM)    as   the   miss   distance    the   missile    cross    the    target    with 

no    commands.  It    can  be   calculated      at    t=0    as      the    initial 

distance    between    target      and    missile    plus    the      miss   distance 

due    to   the    gravity    or: 
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ZEM    =    z    (0)     +    -j"    9   f- 


(1.3*) 


which   is   equal    in  all  the   three   first    cases. 

In  the  first  three  cases,  the  initial  missile's 
commanded  acceleration  is  the  same  as  seen  in  table  3.1. 
Considering  that  the  control  Ac  necessary  to  correct  the 
initial   miss  distance  can  be   caiculited    as: 


Ac 


1 


Zi*< 


Ti 


Ti_ 


=    26.7 


(3.35) 


which   is   close    to   the   initial   control    Ac. 

This  suggests  that  the  initial  Ac,  would  be  that  one 
necessary  to  correct  the  initial  ZEM  in  Z  direction,  which 
agrees  with  reference  1.  Notice  however  that  in  all  cases 
the  initial  Ac  is  less  than  the  calculated  value  of  Aco, 
which  agrees  with  the  previous  statement  that  some  error  was 
expected   in   the    initial   part  of   the   computations. 

Considering  the  ideal  initial  roll  angle  as  defined 
before,    one   has: 


*, 


deal 


=   tan-1 


Z£My 
ZE  M, 


(3. 36) 


From  table  I,  and  figures  showing  roll  angles,  it 
can  be  verified  that  the  missile  is  banking  to  reach  angles 
larger  than  ^(  .  ,  in  order  to  correct  its  trajectory  to 
hit    the   target. 

Therefore, in  all  cases,  the  missile  begins  it- 
trajectory    with    an    Aco    (discussed      before)  ,       and    a    roll   rate 
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which  is  proportional  to  the  targst  acceleration  in  the  ? 
direction.  As  the  missile  rolls  at  decreasing  roll  rates, 
the  commanded  acceleration  is  changed  in  order  for  to 
compensate  the  effect  of  the  roll  angle  on  its  Z  direction 
component. At  some  time  when  the  control  Pc  is  zero  or  near 
zero,  the  control  Ac  begins  to  follow  a  linear  law,  as 
suggested  in  eqn.  2.7: 


Ac=    Aco 


[-*-] 


Notice  however  that  the  term  Aco  in  this  equation  is  no  more 
the  actual  initial  commanded  acceleration,  but  that,  one  the 
missile  would  have  if  its  initial  roll  angle  was  equal  to 
the    final 

Such  behavior  defines  a  boundary  in  the  control  Ac 
which  is  clearly  seen  in  3.29,  where  the  commanded  accelera- 
tion is  bounded  by  the  curve  of  the  control  Ac  the  missile 
would  have  if  its  initial  fia  was  ejual  to  ^o|dea((°rr  if  no 
commanded  roll    was    necessary   to   reach    the   target). 

Althoug  the  missile  commands  roll  angles  larger  than 
the  ideal  initial  roll  angle,  it  is  possible  to  do  a 
prediction-  with  no  no  computer  work-  of  an  aproximation  for 
the  maximum  acceleration  the  missil=  would  experience  during 
its    fligth,    as    folio  wing  (see   fig.    3.32): 


A,,  = 


ZEM2 

z.         a> 

(3.37) 


CO 


s^ 


c/to| 


where    $60         cames   from  eqn.  3.50,      and    Ac      is    a   streight   line 

'      >cital 

as    in    fig   3.32. 

Other  interesting  point  is  that  the  missile  is 
commanding  to  reach  roil  angles  about  twice  of  pe  lde,  when 
the    target      is    at      small   accelerations,         and   when      at    large 


49 


accelerations,  the  missile  is  making  a  small  correction  on 
its    roll   angle,    as    shown   in    table   I. 

One  can  see  from  the  figures  showing  miss  distance, 
that  the  relative  position  of  missile  to  target  is  about  the 
same  in  all  three  initial  cases.  Thus,  the  new  $ldeai  at 
each  point  is  the  same.  ks  the  0O  dcaX  computed  it  t=0  is 
smaller  when  the  target  is  at  snail  accelerations,  the 
correction  has  to  be  larger  in  order  to  hit  the  target,  as 
seen    in    figure    3.30. 

In  table  I,  one  can  see  the  final  miss  distances, 
the  final  miss  distance  eg  to  eg  and  the  final  roll  angle. 
Such  results  show  that  with  this  iigitalized  model,  good 
results   has   been   obtained. 
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Figure    3.1        Representation    of    the   System, 
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Figure   3.2        ?ariation  of  Commands  with   Initial   Roll   Angle. 
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Figure  3.3    Variation  of  Roll  Angle  -  Small  Angles. 

53 


Figure  3.4    Variation  of  Roll  Angle  -  Large  Angles. 
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Figure   3.5        Ideal   Initial    Roll   Angles. 
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Figure  3.7   Commanded  Acceleration-  Case  0. 
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Figure  3-8    Commanded  Roll  Rate-  Case  0. 
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Figure   3.9        Miss  Distance   in  Z    Direction-   Case   0, 
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61 


1ST  CRSE 

:nitjal  target  acceleration-  -.5  5 

initial  target  pos1tion-1qq  m 

sample  per10d-0.05  sec 


n 

1/1 
rs- 

0 
in- 

OJ 

in 
PT 

LJ  ta- 
in R 

O 

IT) 

m  ts- 

CC  - 

UJ 
t— 
U 

3Z  0 

1  £- 

0 

01 
U 

_] 

0  0- 

m 

a 

10" 

m 
cn!" 

o 
Q. 

Q.O   10. 0   20.0   30.0    40.0    SO.O    60.0    70.Q   80.0    90.0   100.0 

K 


Figure  3.12    Commanded  Acceleration-  Case  1 


62 


u  "J 
CO  <= 
\ 

O 

q: 

i  in 

(— 

CE 

a: 


CC.O 


1ST  CR5E 

INIT1RL  TRRGET  ACCELERATION-  -.5  G 

INITIAL  TARGET  POSITIQN-100  M 

SRMPLEPER10D-0.05  SEC 


0.0        10.0        20.0        30.0        40.0 


50.0        60.0        70.0        80.0        90.0        100.0 
K 


Figura    3.13        Commanded   HoLl    Rate-   Cass    1. 
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Figure    3.14        Miss   Distance-   Casel. 


64 


1ST  CASE 

INITIAL  TARGET  ACCELERATION-  -.5  G 

INITIAL  TARGET  POSITION-IOO  M 

SAMPLE  PERI OD-Q. 05  SEC 


o 


o 

in 

a 

mi  ii  ii'inmii  r  I'liimmi 

o 

^" 

a 

at 

o 

CO     •- 

a  ° 
o_ 

i 
u 

_J  m 
CD  ™ 

D 

,     a 

i      D 

!  d 
!  a 
■a 
_i5 

i°' 

!         a ! 

d* 



a 
a 
!     a 
— j— tr— 
!  a 
! D 

4 — 

■[ 

-L 

..L— 

__-  . 

__._-.u- 

1 

in 

;a 
__J3 

a 

a! 

a  ! 

a 
d~ 

1 

C 

a 

3..L 

! 

o 

a 

o 

a*" 

□ 

o 

a 

■q 

0.0        10.3 


20.0         30.0         40.0         50.0         50.0 

K 


70.0         80.0         30.0        100.0 


Figure    3.15        Roll   Angls-   Case    1. 
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Figure  3,16    Commanded  Acceleration-  Case  2. 
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Figure    3.17        Commanded    RoLl    Rate-   Case    2. 
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Figure    3.18        Miss   Distance-   Case   2, 
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Figurs  3.19    Roll  Angls-  Case  2. 
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Figure  3-20    Commanded  Acceleration-  Case  3. 


70 


r— 1 


IP 

a 


...  0. 


cj  m 

CO  t^ 

co 

Q 

cr 

1   °^ 
i—  • 

<r 


o 


a. 


3RD  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  P05ITION-100  M 

SAMPLE  PERI0D-Q.Q5  SEC 


rmrrnn 


BffiWffi 


bp—hm 


0.0        10.0         20.0         30.0         10.0         50.0         60.0         70.0         80.0         90.0        100.0 

K 


Figure  3.21    Commanded  RoLl  Rats-  Case  3, 
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Figure    3.25        Commanded   RoLl    Rate-   Case    4. 
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Figure  3.26    Miss  Distance-  Case  <*. 
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Figure  3.27    Roll  Angle-  Case  4. 
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Figure  3.30    Corrections  on  the  Roll  Angle. 
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IV.    ANALISYS    OF    GAINS,     SAMPLE    RATE    AND    PITCH    ANGLE 

A.       ANALISYS   OF    THE    GAINS 

In  chapter  3  has  been  developed  a  solution  for  the 
optimal   control    of   a    system    as: 

x(k*1)  =   A  (k)     x(k)  +    E(k)    u  (k)  ♦   E   g 

It  would  be  interesting  to  check  if  -he  optimal  gains 
reach  steady  state, but  at  the  momsit  -hat  the  extention  for 
large  roll  excursions  has  been  introduced,  and  the  system  is 
being  feed  with  optimal  commands  which  are  varyina  each  step 
of  time,  such  idea  can  not  be  applyed.  However  one  can  do 
such  check  in  the  model  for  small  roll  excursions,  which  is 
acctually   used    to   compute   the    optimal    commands. 

Doing  this,  one  has  the  time  history  of  those  gains  as 
in    4.1   and    4.2 

One  can  notice  in  fig.  4.1  that  the  gain  FG(2,1)  associ- 
ated with  the  effect  of  gravity  has  no  effect  on  the 
ccmanded  roll  rate,  and  that  FG(1,1),  as  shown  in  fig.  4.2, 
has  a  large  effect  on  the  ooaaanded  acceleration. 
Furthermore,  this  gain  reachs  stsaiy  stats  very  fast,  thus 
it  can  be  assumed  that  the  gain  FG  will  be  egual  tc  the 
steady   state  value    during   all    time. 

From   egn.    3.19,     and   assuming    steady    stats: 


u(k)=    -F(k)    x(k)-    FG(k)    g  (4.1) 


and    substituting   g: 
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u  (k)  =    -F  (k)    x  (k)  -    C 


(4.2) 


where  the  second  term  in  the  right  hand  side  is  a  constant., 
and  its  value  is  exactly  equal  to  tie  value  of  the  commanded 
acceleration  necessary  at  t=0  to  correct  the  gravity  fall  of 
the    missil    (or    to  correct,  the    initial   ZEM   due   to    gravity)  . 

It  might  be  supposed  that  one  could  solve  the  optimal 
control  problem  for  the  system  represented  by  eqn.  4.1,  just 
considering   one    reduced    system   represented    by: 


x(k+1)  =   A  (k)     x(k)  +    B(k)    u  (k) 


(4.3) 


with    a    bias   in    the   ccntrol    as: 


u  (k)=   -F  (k)    x  (k)  ♦    C 


(4.4) 


But      as      shewed    in      the      following      analysis,      this      is      not 
possible. 

The    constant    term  in   the   right   hand    side    of    eqn.      4.4    is 
calculated    as   follows: 

from   eqn.: 

initial    ZEM    due    to  gravity=   ZEM   =    --    g   t2 


C    = 


£M 


r   *2 


1  t      I  J 

---     -      ---  CCS   0 

.     I  6Ti   J  ~° 


and    the    gains   F(k)       are      calculated    using   a    Ricatti    equation 
as    usual. 
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Case  5  was  tested  using  the  above  considerations,  and 
using   the   same    scenario   of    case   3. 

Figures4.3  and4.4,  show  the  time  history  of  the 
controls.  One  can  see  that  the  commanded  acceleration  has 
begun  at  same  values  as  in  case  3,  but  the  commanded  accel- 
eration reachs  a  peak  considerable  higher,  them  dcreases 
does  not  following  a  linear  law,  with  a  final  at  14.7  meters 
per  second  sjuare,  being  this  terminal  value  due  to  the 
constant   term  representing    the   effect    of   the    gravity. 

Refering  to  fig.  4.4,  the  commanded  roll  rate  begins  at 
a  same  value  as  in  case  3,  but  ae  the  control  Ac  is  too 
high,  it  reachs  negative  values,  going  to  zero  almost  at  the 
end  of  the  running  time.  This  behaviour  of  the  control  leads 
to   a    large    miss    distance   as    seen    in   figure    4.5,    and   tablell. 

Figure  4.6  shows  the  time  history  of  the  roll  angle, 
which  rises  to  values  close  to  1.5  radians.  As  the  accelera- 
tion at  this  point  is  larger  than  the  correct  value,  the  the 
corrections  are  excessive  and  the  roll  angledecreases  at  the 
end    of   the    running    time    to    the    value    of    .12   radians. 

Thus,  one  can  see  that  the  gain  due  to  the  gravity's 
acceleration  can  net  be  replacai  by  its  steady  state 
value. This  kind  of  simplification  can  thus  not  be  done  in 
the    system    being  studied. 

B.       EFFECT    OF    THE   SAMPLE    RATE 

Throughout  all  the  simulatios  a  sample  period  of  .05 
seconds  has  been  used. In  this  section  a  brief  study  of  the 
effect   of   the   change    of   this   sample    rate    is   giv=n. 

Two  best  cases  have  been  selected  to  ilustrate  the 
effect  of   the   sample    period. 

The  first  case,  case  6,  has  been  run  with  a  sample 
period  of  .1  seconds  and  consits  Df  the  same  scenario  as 
case    3. 
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As  one  can  see  in  figure4.7  and  tablelll,  there  is  no 
noticable  change  in  the  commanded  acceleration,  but  the 
commanded  roll  rate  begins  at  smaller  value  than  in  case  3r 
as  is  seen  in  fig. 4. 8  This  initial  decrease  in  roll  rate, 
leads  to  a  large  miss  distance  in  Y  direction  as  shown  in 
fig. 4.9,    and  to    a  small   value    of    roll    angle    (se  =    fig. 4.  10). 

The  second  case,  case  7, has  a  smple  period  cf  0.025 
seconds.  There  is  no  noticeable  change  in  the  time  history 
of  the  control  Ac  as  shown  in  fig. 4.  11.  The  commanded  roll 
rate  begins  at  a  higher  value  than  in  case  3  as  shows  fig. 
4.12,  which  leads  to  a  final  miss  distance  in  Y  direction  of 
-2.22  meters  and  in  Z  direction  smaller  than  case  6  (see 
fig4.13).  Figure4.14  shows  that  the  final  roll  angle  is 
increased  and  the  missile  cross  the  target  wi-h  1.28  radians 
and    with   a    CG-to-CG    distance  of   2.5    me.ers. 

In  both  cases,  the  miss  distance  was  increased  over  the 
nominal  value  obtained  with  a  sampling  rate  of  0.05  seconds. 
Thus,  it  would  appear  that  there  is  an  optimal  value  for  the 
sampling  rate,  which  may  be  conected  wi-h  the  geometry  of 
the    scenario  and   with  time    to   go. 

C.       EFFECT    OF   THE    INITIAL    PITCH    ANGLE 

It  is  important  at  this  step  to  remenber  that  througout 
this  work  has  we  have  been  discussing  a  dimensional  model, 
where  there  is  no  information  on  the  X  coordinate,  so  it  is 
impossible    to   verify    the   behavior   of    the    pitch   angle. 

In  this  work,  since  in  all  the  previous  scenarios  the 
initial  angle  0  was  equal  to  zero,  this  value  has  been  kept 
as  a  constant  during  all  time,  and  considering  that  whithout 
any  information  of  a  third  dimension  it  was  not  possible  to 
correct  the  time  to  intercept,  this  time  was  also  kept 
constant    and  equal    tc  the    nominal    value    of    5    seconds. 
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Notice  -hat  this  assumption  is  likely  to  be  correct  if 
one  has  the  horizontal  initial  iistance  from  target  to 
missile  compared  with  the  initial  vertical  distance  between 
target  and  missile  large  enough  in  order  tc  have 
small        angles. 

The  guestion  that  rises  is,  how  could  this  pitch  angle 
affect  the    system  if   it    was    net   small? 

As  seen  in  fig. 4. 15,  the  missile  velocity  in  the  X 
direction   would    be: 


Vmx  =   Vm  cos©   ♦   g    cos6>  sin©    t  (4.5) 


which  has  an  affect  not  only  from  tha  pitch  angle,  repre- 
sented by  the  cos  0  ,  but  also  from  tha  gravity's  accelera- 
tion,   which   will   leads  to   a    different    time   to   intercept. 

Considering  the  same  physical  soenario  as  in  case  4,  but 
changing  the  initial  pitch  angle,  in  order  to  have  the 
missile  pointing  to  the  target  (see  fig. 4. 16),  and  keeping 
tha  missiles  velocity  of  10  90  m/sec  in  the  X  direction,  one 
has  a  ccmpletly  different  geometry  of  the  problem  as  seen 
from    tha   flight    path    reference   frame. 

With  this  new  situation  (see  fig.  4.  16)  ,  case  8  has  been 
run.  Figures  4.17  and  4.18  show  that  now  the  missile  is 
comanding  large  positive  accelerations,  and  the  roll  rate  at 
the  begining  of  the  flight  is  too  high,  going  to  zero  in  a 
very  small  period  of  time.  The  miss  distance  as  seen  in 
fig. 4.19  are  increased  in  the  initial  part  of  the  fligth  and 
as  the  missile  corrects  its  trejectory  it  is  decreased  to 
reach  a  final  CG-to-CG  miss  distance  of  2.17  m.  The  roll 
angle,  due  to  the  large  control  Pc  is  oscilatory  in  the 
begining  and  becames  constant  with  a  value  of  .57  radians 
(see    fig. 4.  20   and   table    IV)  . 
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Notice  -hat  the  high  values  of  acceleration  needed  are 
in  some  part  dae  to  a  vertical  component  of  target's 
velocity,  which  is  seen  from  the  flight  path  frame  as  the 
target  was  raanouvering  in  the  Z  direction  with  constant 
velocity.  These  large  accelerations  leads  to  roll  rates  too 
large  for  the  physical  integrity  of  the  missile.  This  means 
that  although  althoug  the  good  results  obtained,  if  compared 
with    case   4,   they   are  not   practical. 

In  order  to  get  rid  of  ths  vertical  manouever  of  the 
target,  case  9  has  been  run.  In  this  oase,  the  scenario  is 
the  same  as  before  with  the  target  also  pointing  down,  with 
the  same  pitch  angle  as  the  missile,  and  has  a  X  velocity 
equal   to   the  previous  case    (see   fig. '4. 16). 

From  fig. 4. 21,  one  can  see  that  the  decrease  in  the 
control  Ac  is  substantial  if  compared  with  case  8,  but  the 
commanded  roll  rate  is  still  toe  large  as  shows  fig.  4.22  Ths 
time  history  of  the  miss  iistancs  is  better,  resulting  in  a 
final  CG-to-CG  distance  of  about  50^  of  case  8  (see 
fig. 4.23).  The  roll  angle  is  not  oscilatory  as  seen  in 
fig. 4.24  and  the  missile  oross  the  target  with  1.3  radians 
in   roll. (see  tablelV) . 

The  results  obtained  in  the  two  previous  cases,  suggests 
that    the      algorithm    developed    in      this    wor!c    could      be   readly 


applied    to      air-to-surface    missiles 


:n    the      latter, 


:he 


scenario  would  be  favorable  to  the  missile  than  in  either 
previous  cases,  since  one  can  consider  that  the  target  could 
be  essencially  stationary  in  comparison  with  -.he  missile 
speed. 

Case  10  has  been  run  with  this  assumption,  and  the  scen- 
ario as  in  fig.  4.16.  The  target"  is  with  zero  velocity  and 
acceleration,  and  the  missile  begins  its  flight  600  meters 
above    it,    with    ths    same    initial   pitoh    angle    as   before. 

Fig. 4. 25  shows  the  time  history  of  the  commanded  accel- 
eration,     where   one    can   see    that   as   there   is    no    roll   rate   to 
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be  commanded,  the  acceleration  is  following  a  straight  line 
with  very  resoanable  values.  The  miss  distance  is  show  in 
figure4.27,  which  shows  the  final  3G-to-CG  distance  of  .31 
meters. 

Based  in  the  results  of  these  tests,  one  can  see  -chat 
there  will  be  some  effect  of  the  pitch  angle  on  the  miss 
distance,  not  only  due  to  its  effect  on  the  time  to  inter- 
cept, but  also  because  at  the  moment  that  there  is  a  pitch 
angle  different  from  zero,  even  if  the  Target  is  keeping  its 
flight  level,  in  the  flight  path  frame  a  component  of  the 
target's  velocity  will  show  up  leading  the  missile  to 
command   large      accelerations   and    roll    rates.  Althoug   this 

results  harm  the  performance  in  an  air-to-air  missile,  in 
the  case  of  air-to-surface  missiles,  when  the  target  has 
been  considered  with  no  motion,  good  results  has  been 
obtained. 

D-       EFFECT    OF   TIME    TO   INTERCEPT 

In  the  simulation  of  case  1  aid  2  in  chaptar  3, it  has 
been  observed  that  when  the  target  was  at  small  accelera- 
tions, the  missile  did  larger  corrections  on  its  roll  angle, 
with  respect  to  the  ideal  initial  roll  angle,  than  when  the 
target  was  with  large  accelerations.  One  can  think  that  must 
be  seme  kind  of  compromise  between  the  velocity  rate  of 
target  and  missile  (which  will  reflected  on  the  time  to 
intercept),  and  the  relative  position  between  them,  which 
will    affect   the    miss    distance. 

In  order  to  do  a  brief  analisys  on  this,  case  11  and  12 
has    been   run. 

In  case  11,  the  scenario  of  case  2  has  been  kept,  with 
the  exception  that  the  missile's  velocity  was  change  to  2000 
m/sec,    which   means    that    Ti    was   changed    to    2.5    seconds. 
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Figure  4.29  shows  that  the  icceleration  is  largely 
increased  due  to  the  small  time  required  to  correct  the  ZEM, 
and  the  commanded  rcll  rate  is  almost  twice  of  case  2  (see 
fig.  4.30).  The  final  miss  distance  is  more  than  four  times 
the  value  obtained  in  case  2,  as  seen  in  fig.  4.31  and  table 
V.  The  final  roll  angle  is  about  half  as  in  case  2,  since 
the  projected  final  position  of  the  target  in  the  Y  direc- 
tion   is   less  than  in   case  2     (see    fig.    4.32). 

Case  12  was  run  with  a  scenario  less  favorable  to  the 
missile,  where  all  the  conditions  of  the  previous  case  was 
kept,  except  the  target  position  that  has  been  increased  to 
200    meters    above   the    missile. 

Now,  the  commanded  acceleration  are  much  larger,  with  a 
initial  Aco  of  103  m/sec2,  being  almost,  impossible  to  see 
the  difference  of  the  time  history  of  the  acceleration  from 
one  straight  line,  as  shown  in  fig.  4.33  The  commanded  roll 
rate  is  small,  about  the  same  as  in  case  2  (see  fig.  4.34). 
The  change  in  the  miss  distance  is  noticeable,  with  a  final 
CG-to-CG  distance  of  4.8  meters  as  in  table  7  ,  and  figure 
4.35.  The  final  roll  angle  explain  the  shape  of  the  acceler- 
ation curve,  since  with  the  small  roll  angle  as  shown  in 
fig.  4.36,  the  system  is  behaving  as  for  small  roll  excur- 
sions. 

Notice  that  from  this  analyzis,  one  has  to  realize  that 
there  is  seme  kind  of  envelope  where  the  2-D  system  is 
valid.  And  in  order  tc  determine  this  envelope,  one  has  to 
take  in  account  not  only  the  geometry  of  the  scenario,  but 
also  the  time  to  intercept,  which  is  determined  not  only  by 
the  relation  of  velocities  of  missile  and  target,  but  the 
pitch   angle   too. 
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Figure   <4.  1        Gains   affecting   the   commanded   acceleration 
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Figure   4.2        Gains  Affecting   the   Commanded   Roll   Rate, 
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5TH  CASE 

INITIRL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  POSITIQN-IOQ  M 

SAMPLE  PEPJOQ-0.05  SEC 


O 

... J&Ji..      '  ... L ! L. 

o 

o 

8" 

o 

*\  \\               i        ! 
r % r r r 

a                a 

1             o 
[            •  a 

i            la!            !            ! 

a     :         !      a :         : 

o 

1     1     fc    i     i 

.JS. L i. LB»— 1. L 

!            !     a      .            . 

i        Q     i                • 

J            !            !        a  !            : 

XL  o 

H              i                  i                  .                 CD                 • 

!                          P           ! 

,                                        o 

1    £" 

2° 

p   " !        1        1        :  \  1 

£8" 

1, 1 

LJ  tn- 
O  N 
o 
cc 

o 

J           !            !            !            !            !     tj. 

._._  JH3k«. ! ___L_____. 

8" 

o 

o 

o 

o 

i      '  i        l  '   '  i  '      i — H 

f         i    »'    1      i[ 

0.0        10.0         20.0         30.0        10.0 


so. a 
K 


50.0         70.0        00.0         90.0        100.0 


Figars    4.3        Commanded   Accaleration-Case    5, 
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5TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 
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Figure   <*.4        Commanded    Roll    Hate-Case   5 
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5TH  CRSE 

INITIBL  TRRGET  RCCELERPTION-  -4.  G 

INITIAL  TARGET  POSITION-lOO  M 

SfitfPLE  PER1QD-0.Q5  SEC 


Figure    4.5        Hiss   Distanca-Case   5. 
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5TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  POSITIQN-lOO  M 

SAMPLE  PERIDD-O.Q5  SEC 
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Figure    4.6        Roll   Angle-Case   5- 
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TABLE   II 
Results   from   Biased   Control 
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6TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  P05ITI0N-100  M 
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Figurs    4.7        Commanded   Accslsration-Case   6. 
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6TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  POSITION-1QO  M 
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Figure   4.8        Commanded   Roll    Rate-Case   6. 

93 


6TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  PCSITION-lOO  M 

SAMPLE  PERIOD-O.10  SEC 
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Figure    U.9        Miss   Distaaca-Case   6, 
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6TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  P05ITI0N-100  M 

SAMPLE  PERI OO-O. 05  SEC 
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Figure    4.10        Roll    Angls-Case   6. 
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7TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  P051TI0N-1Q0  M 

SAMPLE  PERIQD-0.025  SEC 
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Figurs    4.12        Commanded    Soil   Rate-Case   7 
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7TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  POSITION-100  M 

SAMPLE  PERI0D-0.Q25  SEC 


Figure    4.13        Miss   distiaca-Case   7. 
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7TH  CASE 

INITIRL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  PQ5ITI0N-10Q  M 

SAMPLE  PERIOD-0.025  SEC 
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Figure    4.14        Roll    Angle-Case   7. 
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Figurs    4.15        Effect   cf   Initial    Pitch   Angle. 
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Figura    4.16        Scenarios    With    Pitch   Angls 
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8TH  CRSE 

INITIAL  TRRGET  ACCELERATION-  -4.   G 

INITIAL  TRRGET  POSITION— 600  M 

SAMPLE  PERIQD-0.05  SEC 


Figurs    4.17        Commandad   Accaleration-Case    8. 
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Figure    4. 1 8        Commanded    RoLL    Rate-Case    8. 
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8TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  POSITION— 600  M 

SAMPLE  PERIOD-0.05  SEC 
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Figure    4.19        Miss   Dis-canca-Case   8. 
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8TH  CR5E 

INITIRL  TRRGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  POSITION— 500  fl 

SAMPLE  PERlOD-0.05  SEC 
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Figure    4.20        Roll    Angla-Case   8. 
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9TH  CP5E 

INITIAL  TARGET  ACCELERATION-  -4.   G 

INITIAL  TARGET  POSITION— 600  M 

SAMPLE  PERIQD-0.Q5  SEC 


Figurs    4.21        Commanded   Accaleration-Cass    9- 
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9TH  CR5E 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  POSITION— 600  M 

5AMPLE  PERIOD-0.05  SEC 


Lfl- 

o 

^"» 

a 
pi- 
es 

o 

o 
d- 

y>  » 

\n 

82- 

»-  o 

<n    •« 

CK.  •"< 

_J 

-1  o 
O     ■. 

sc  * 

o 

oi" 

o 

o 
in" 

o 
ri" 

o 

1-%" 

d. 

i                                   j            j                        j            ;            ;                        ; 
i            i            i            i            i            i            i            i            i            i 

a       i        ;        ;        :               ;               :        :        : 
*        i        *        i        *        i        i        i        i        t 

!          !          *          !          !          I          !          f          I          I 
i          i          i          »          i          i          i          i          i          i 

!!*!!!!!!! 

i          i          i          >          i          <          i          i                    i 

i          i          i          >          i          i          i          i          i          i 
i          .          i          •          .          .          i          i          i          i 

!                      !           I           I           I           '           !           I           I 

a        ;           ;           ;           ;           ;           ;                      ;                      ; 

d      ;           ;           ;           ;                                           ;                     j 

0.0       10.0       20.Q       30.  □       10.0       50.0       60.0       70.0       80. Q       90.0       100.0      110.0 

K 


Figure    4.22        Commanded   Roil    Rate-Case   9 
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9TH  CR5E 

INITIAL  TRRGET  ACCELERATION-  -4.  G 

INITIRL  TRRGET  POSITION— 500  M 

SPMPLE  PERI0D-Q.C5  SEC 
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Figure    4.23        Hiss   Distaaca-Case   9, 
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9TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  POSITION— 500  M 

SAMPLE  PER10D-0.05  SEC 
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Figure    4.24        Roll    Angle-Case   9. 
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10TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 

INITIAL  TARGET  POSITION— 500  M 

SAMPLE  PERI0D-0.Q5  SEC 
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Figure  4.25    Commanded  Accslaration-Case  10. 
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10TH  CASE 

INITIAL  TBRGET  ACCELERATION-  -4.  G 

INITIRL  TRRGET  POSITION— 600  fl 

SAMPLE  PERJQD-Q.05  SEC 
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Figure    4.26        Commanded   Roll   Rate-Case    10, 
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!07H  CASE 

INITIRL  TRRGET  ACCELERATION-  -4.   G 

INITIAL  TARGET  POSITION— 600  M 

SAMPLE  PER  1 00-0. 05  SEC 


a 

a 

o 

K  o 

U 

1 
LJ 

^ 

!a 
a 

1°' 
IT) 

Q 

lrf< 

.0       2^.Q       2}. 

a     39.0     3! 

k    : 

1 
i 

3aaoODnnat 

1 
i 
i 
1 
1 
t 
1 

1.0       1<J,0       1! 
rfrtjD?° i 

■;S°°* 

1,0      5 

E 

a 

7°*      . 

i 
i 
1 

■ 
1 

■ 

i 
o 

iiiii 

iiiii 
1      1      1      1      1 
iiiii 
1      1      1      1      1 
1      1      1      1      1 

1 

1           >           1 
1           1           1 

O 

LEGEND 
a  -  MISS  DIST.   Z  DIRECTION 
0-  MISS  DIST.   Y  DIRECTION 

u 

* 

O 

6 

1                       1                       1                       1                       c 

1                 1                 1                  1                 1 

Figure    4.27        Miss   Distanca-Case    10, 
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10TM  CR5E 

INITIBL  TARGET  ACCELERATION-  i.   G 

INITIRL  TRRGET  POSITION— BOO  M 

SRKPLE  PERIQD-0.Q5  SEC 
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Figure    a. 28        Roll    AngLs-Case    10. 
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nsi  crsc 

INITIflL  TPPGET  RCCELEPBTION-  -1.  Q 

INITIP.L  TRRGE7  POSITION- ICO  M 

SRMPLL  PER iOO-O. 05  SEC 
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Figure  4.29    Commanded  Acceleration-Case  11. 
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11  ST  CH5E 

INITIAL  rfiRGET  RCCELERRTIGN-  -J.  C» 

INITlfiL  TPRSST  POSITION- 1QC  M 

SRMPLE  Pl^IOC-O.GS  5lC 
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Figure    4.30        Commanded   Roll   Rate-Case    11 

122 


!15T  zr^L 

INITIAL  WSET  FCCELDWTJON-  -1.  G 

INITIAL  TARGET  POSITION- 100  M 

SFtfPLE  PER)GD-C05  5CC 
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Figure    4.31        Miss   Distaaca-Case    11. 
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U5T  CPSC 
INITIAL  TBR5F:T  ACCELCRATION-  -1. . 
INITIAL  TBRGST  POSITION- 100  K 
SAMPLE  PER10D-0.Q5  SEC 
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Figure    a. 32        Roll   AngLs-Case    11. 
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12TH  CASE 

INITIAL  TARGET  ACCELERATION-  ~i.  8 

INITIAL  TARGET  PQSITION-20Q  M 

SAMPLE  PER10D-O.Q5  5EC 


o 


O-  .- 


UJ  o-  • 

o 

tr.  ^ 

u 
u 

.  8' 


Ui 


a 

-■a- 


r-tJ- 


I 1- 


.lD.. 


a— ■ 
a 


■rpp- 

a 


I     I     I 1 1 1     I     I     I 

0.0    S.O    10.0    15.0    20.0    25.0    30.0    35.0    10.0    15.0    50.0 

K 


Figure   4.33        Commanded   Accaleration-Case    12. 
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12TH  CASE 

INITIAL  TARGET  ACCELERATION-  -1.  G 

INITIAL  TARGET  P0SITI0N-20Q  M 

SAMPLE  PERIOD-0.05  SEC 
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Figure  4.34    Conianded  fioLl  Rate-Case  12. 
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12TH  CRSE 

INITIRL  TARGET  RCCELERRTION-  -1.  G 

INITIBL  TRRGET  P0SITI0N-2D0  M 

SAMPLE  PER1QD-0.Q5  SEC 
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Figure    4.35        Miss   Distancs-Case    12. 
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12TH  CASE 

INITIAL  TARGET  ACCELERATION-  -1.  G 

INITIAL  TARGET  P0SITIQN-2QQ  II 

SAMPLE  PERIOD-0.05  SEC 
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Figure    4.36        Eoll    Angls-Case    12 
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?-     ZINAL   CONCLUSIONS    AND    COMMENTS 

The  scops  of  the  present  work  was  the  development  of  an 
optimal  digital  control  to  be  applied  on  a  bank-to  turn 
missile. 

A  two  dimensional  model,  as  sugssted  in  reference  1,  was 
adopted.  After  the  digitalizat ion  of  the  continuous  model  it 
was  necessary  to  solve  a  modified  Ricatti  equation  since  in 
the  stats  equation  there  was  a  third  term  representing  the 
gravity's  effect.  The  approach  that  has  been  adopted  is  new, 
and  although  good  results  were  obtained  for  the  scenarios 
considered  in  this  work,  is  necessary  that  the  algorithm  be 
further  tested  and  evaluated  in  similar  problems  due  to  its 
novelty. 

The  optimal  was  solved  with  an  initial  restriction  to 
small  angles.  This  condition  was  later  relaxed  so  that  large 
roll    angles   could   be   analyzed. 

It  is  dificult  to  compare  the  present  work  with  previous 
results  since  Stallard  has  indicate!  a  mistake  in  his  orig- 
inal   paper,    and    further   works    in    this    area   was   net    found. 

However  some  comparison  with  Stallard  work  is  possible. 
Tha  commanded  acceleration  cf  the  missile  are  such  as  to 
correct  the  ZEM  at  each  point,  this  agrees  with  that  refer- 
ence.There  is  a  proportional  relationship  between  the 
commanded  roll  rate  and  the  commandad  acceleration  ,  and  the 
commanded  rcll  rate  is  proportional  to  the  defined  $ ldea/  at 
each    point,    which   again    agrees    with    reference    1. 

The  algorithm  develops!  in  this  work  requires  extensive 
computation  at  each  step,  and  is  olear  that  some  software 
optimization  will  be  needed. The  motivation  for  considering 
the  constant  steady-state  gain  due  to  gravity  was  to 
decrease  this  computational  burden.  This  approach  however 
resulted  in   unacceptable   miss   distance. 
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Another  point  of  investigation  that  could  reflect  or.  the 
period  available  to  the  computer  to  perform  its  calculations 
was  a  change  in  the  sample  rate.  Itfo  different  sample  rates 
were  investigated,  both  lead  to  larger  errors  than  the 
nominal  period  of  .05  seconds.  A  detailed  study  on  this 
issue  is  'left  as  sugestion  for  future  works,  since  some 
optimal   value   of   the   sampling   rate   is    clearly    indicated. 

It  is  important  to  keep  in  mind  that  the  model  adopted 
is  two  dimentir.al,  while  the  actual  problem  is  three  dimen- 
sional, thus  some  brief  studies  ware  conducted  in  order  to 
check   the   region   of    validity   of   tha    2-D. 

In  the  analysis  of  the  pitch  angle,  one  car.  see  that  is 
necessary  to  have  small  variations  in  pitch  in  order  to 
approximate  it  as  a  constant.  However,  at  the  moment  that 
this  angle  is  different  from  zero,  as  explained  in  chapter 
4,  it  is  possible  to  have  in  the  flight  path  reference  frame 
a  target  manouver  in  the  Z  direction  that  will  lead  to  large 
acceleration  oommands,  leading  the  missile  to  larae  miss 
distances,  when  considering  a  movable  target.  When  the 
present  system  was  tested  tested  against  fixed  targets,  the 
results  were  guite  gcod,  this  suggests  the  application  of 
the    model   in   air -to-surf ace    missiles. 

Further  investigation  were  mads  on  the  effect  of  time- 
to-go.  As  expected,  decreased  time  to  go,  results  in 
increased      miss      distance.  A   de-ailed      analysis      of      more 

complex  scenarios  is  needed  in  order  to  properly  define  the 
effect   of   time    to  go. 

Also,  it  would  be  interesting  to  extend  the  model  to 
three  dimensions  and  include  the  effects  of  lags  on  the 
system  in  future  works.  Finally,  in  appendix  A,  the 
computer  model  used  in  this  work  is  enclosed.  Some  improve- 
ments in  this  program  can  be  done,  mainly  in  the  data  intro- 
duction,   and  in    some   optimization   of    the   running    time. 
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&PPENDII   k 
FORTRAN    PROGRAM 

These  appendix  provides  a  listing  of  the  computer 
program    used  in    the    present    study. 

Since  the  routines  used  are  non-IMSL,  and  a  small  change 
to  double-pr=cision  was  necessary,  they  are  also  being 
provided. 
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